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Abstract: 

Correlations of two like charged pions emitted from a hot and charged spherically expanding 
nuclear system are investigated. The motion of the pions is described quantum mechanically 
using the Klein-Gordon equation which includes Coulomb field and pion absorption. Flow 
modifies the radial distribution of the source function and rescales the pion wave functions. 
The radii extracted from the correlation functions are calculated in sideward and outward 
direction as a function of the pair momentum. Comparison with recent measurements at 
SIS and AGS energies is made. 
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I. INTRODUCTION 


Measurements of correlations of pions and kaons as a function of their relative momenta 
are used to determine the source radius in heavy ion reactions at intermediate and relativistic 
energies. On the basis of the Hanbury-Brown and Twiss effect(see the correlation 

function is related to the source radius via i?o = hc/Aq, where A q is the width of the 
correlation function. However it is well known that the correlation is a complicated 
function not only of the relative momentum but also of the total momentum of the pion 
pair. The correlation is determined by the size of the region from which pions are emitted 
with roughly the same momenta. This has the consequence that for collectively streaming 
matter this region is smaller than the total source due to the strong correlation between the 
momenta and the emission points of the particles. A further modification of the apparent 
source size arises if pions are often rescattered within the source. Thus, for pions the source 
is opaque |^|, and their emission points lie within a thin surface layer. 

On the other hand the central Coulomb force changes the momenta of the particles 
while moving towards the detector. This latter effect was investigated in refs. |^.^|. Pions 
with small momenta are mostly influenced by the Coulomb forces which act quite differently 
on positively and negatively charged pions. Indeed quite different radii have been observed 
recently in collisions of Au on Au at bombarding energies of 1 GeV || and 11 GeV per 
nucleon || at the SIS and AGS accelerators, respectively. 

The aim of this work is to study the combined effects of central charge, opaqueness 
and flow on the extracted radii. As a model we consider the emission of pions from a hot 
spherical nucleonic system the expansion of which can described by radial hydrodynamical 
flow. To ease the necessary numerical calculations we use a spherical spatial distribution 
while for the momentum distribution a relativistic Boltzmann distribution is used. Starting 
from a covariant formalism for two particle emission (Sect. II) we derive the source function 
in Sect. Ill and obtain a concise expression for the matrix element for the emission of two 
particles. Once we have obtained the matrix element the standard technique is applied 
to calculate the correlation function. Coulomb field and opaqueness are included via mean 
fields in calculating the distorted waves for the pions. Numerical results are finally discussed 
in Sect. VI. 


II. BASIC EQUATIONS 

Here, we briefly review the formulation of the Hanbury-Brown and Twiss effect which 
allows to incorporate the mean fields between the source and the two emitted mesons. At 
asymptotic distance the two mesons move with momenta p and p'. Each meson is described 
by a wave function ip p {x) which satisfies an equation of motion which contains the mean 
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field, and p = (u>, p) and x = (f, r) denote |^] the four-momentum and the position in time 
and space, respectively. It is customary to introduce a source term J(x) from which the 
wave function if) p is generated. With these definitions one can express fi]]To) the probability 
for the emission of two meson in terms of the source operators J + (:c), J{x ) as 

cuu / dN ~ [ d, A x\d A X2d A x 3 d A x 4 x 

dpdp' J 

ljj*(x 1 )'llj*,(x2)lpp’(x3)ll>p(x i )(J + (x 3 )J + (X4 : )J(x 1 )J(X2)). (1) 

The basic assumption for applying interferometry to nuclei is the chaoticity of the 
source i.e. the absence of initial correlations between the two emitted pions except those 
correlations coming from the Bose-Einstein statistics. Thus, one writes 


{J + (x 3 )J + (xi)J(xi)J(x 2 )) = 

(J + (x A )J(x 1 )){J + (x 3 )J{x 2 )) + (J + {x 3 )J{x 1 )} {J + {x 4 )J(x 2 )) ■ (2) 


This allows to express Eq.(d) as products of density matrices 

~ p{p,p)p(p',p')+ I p(p,p') I 2 (3) 

with 


p(p,p') = J j d A xd A x' il’*(x)4>p'(x , )(J + (x')J(x)) 


( 4 ) 


which contains the one-body source term. This term can be connected to a classical source 
function S(p,x) via a Wigner transformation 

<J+(/)J(*)) = J . (5) 

The source function S(k,x) describes the creation of a meson at space-time x with four 
momentum k. 

The problem is further essentially reduced by assuming that the interactions do not 
depend on time. This assumption allows the use of stationary solutions for ip p (x) = 
exp(—iujt) V’p(r) with ui = p° = m 2 + p 2 . Inserting the definition(|j) into Eq.([|) and 
integrating over the variable t — t' we obtain 


P{P,P') = J J drdr’ (r)V’p'(r') 

X / (^3 elk{r_r,) J dt e s[{ 


+ n r + r ' u 

—-—,k), («, —-—)]. 


( 6 ) 


In the interference term for p ^ p' the energy of the pion in the source function is fixed to 
the mean energy of the observed pions. In the case of a thermal momentum distribution 


* The convention % = c = k B oitzmann = 1 is used. 
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without any correlation between space and momentum the integration over k introduces 
correlations between r and r' in the order of the thermal wave length \j2-n / (mT). If plane 
waves are used the integration over the difference r — r' can be carried out which fixes the 
k momentum to half the pair momentum leading to the standard form 0 of the matrix 
element. 

Now, the correlation function C 2 is defined as the two-particle emission function nor¬ 
malized to the product of the one-particle emission functions which is given by the first 
part in Eq.(|^). It is convenient to introduce the average pair momentum K and the relative 
momentum q via 

K = i(p + p'), q = p-p' (7) 


Then, the correlation function reads 


C 2 ( K,q) 



\p(K ~~ 2’K + f)| 2 

lK-*)p(K+lK+i) ■ 


( 8 ) 


For completeness we mention that the final state interaction |ll|,|l2[ between the two 
outgoing mesons has not been considered. Inclusion of this interaction would need replacing 
the product of the two outgoing waves (x') in Eq.(,lj) by a correlated wave function 

\k( x,x'). Including the interaction with the source the function ty(x,x') is, however, the 
complicated solution of a genuine three body problem. As a first approximation one could 
assume that the source interacts mainly with the center-of-mass of the pair while the final 
state interaction is a function of the relative distance between the outgoing particles only. 
As a consequence the function T could be split up into a product leading to the standard 
expression 


C 2 (K,q)f instate = P{q) C 2 (K,q) 


(9) 


where the factor P{q ) is in the simplest case the Coulomb penetrability. In the case of 
Coulomb interaction such an approximation means that the quadrupole and higher momenta 
which influence the relative motion of the pair are neglected. Alternatively, one can try to 
factorize the wave function T as it has been done in calculating proton-proton correlations 


III. HYDRODYNAMICAL PICTURE 

Now we restrict our selves to a simple situation. We consider a hot spherical source from 
which pions are emitted. Since the source expands the problem is not stationary. However, 
we circumvent this difficulty by considering only mesons which move with sufficiently large 
velocities. Those mesons are essentially outside the source and always feel a time independent 
Coulomb potential. Having in mind a collision of Au on Au nuclei at a few A-GeV generating 
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a source of a temperature of about 100 MeV, the thermal velocity of the protons in the source 
is about 0.4 c embedded in a flow field of an average velocity of about 0.3 c. Thus, the model 
may be applicable for pion velocities larger than 0.7 c corresponding to pion momenta above 
100 MeV/c or kaon momenta above 300 MeV/c. However, if the particles are released in the 
center of the source they stay a while within the decreasing part of the Coulomb potential. 
Within our model the only way to treat this problem is diminishing the central charge in 
a heuristic manner. Such modifications of the effective charge have been observed @ in 
describing the 7r“/7r + ratio. 

For ultrarelativistic energies the system does not expand spherically but is preferen¬ 
tially stretched in longitudinal direction. Since our treatment neglects this dimension we do 
not consider the longitudinal correlation function. Due to the longitudinal expansion the 
Coulomb force decreases with time. As an approximate measure of the Coulomb action one 
can replace the central charge Z with twice the rapidity density 2(dN + /dy — dN~ jdy) of 
the net charge as was shown in ref. El- 

In the hydrodynamical approach the momentum distribution is defined by a local tem¬ 
perature T and a velocity field given by the four-vector Here we use typical parameters 
for a heavy ion reaction: constant temperature T, radial mean velocity (/?) and radial size 
Rq. Thus, we start with the source function 
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( 10 ) 


The thermal distribution of the momenta within a fluid cell is coupled via ku = k°u° — ku 
to the four-velocity of the cell. For large pion density one should replace the Juttner 
or relativistic Boltzmann distribution with a Bose distribution. A spherically expanding 
system can described by the flow velocity field 


= (it 0 , u), u° = y/l + u 2 , u = f3 0 r/R 0l 


( 11 ) 


where the four-velocity scales with the distance from the center. The parameter /3q is 
related to the mean flow velocity (/3) = (|u|/u°) averaged over the density (© which leads 
to (P) = y/8/nPo in the limit of small /3q. Since fio characterizes a four-velocity it does not 
have an upper bound while in contrary (/?) never exceeds the velocity of light. In Eq. @ 
the pions are radiated off during the emission time r. 

Inserting the source function (0) into Eq.(jiij) and integrating over the time we obtain 

p(p,p') = ( 27r ) 9 / 2 fl 3 e~ (a) ~ a) ' )2 ' r2/2 J drdr'i/jp(r)ip pl {r') 

x J dke _ik ^ r_r, - ) e kU//T e _ ^ r+r, - )2 ^ 8 ' R °^ ( - a ' +w, ^ 1+u2 ^ 2T ^ , (12) 

where the velocity field u is proportional to i -(r L . In the range of the convergence of this 
integral the integration over k leads to a 6 function S(r — r' — iu/T) which allows carrying 
out the integration over (r — r'). Thus, the final expression reads 
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1 


p(p,p') 


(27T 


_l_ e (« u ,') 2 t 2 /2 J dr[i/j p {ar)}*ip p ,(ar) 


xe 


-r 2 /(2flg)-(^+w')V 1 +(/ 3 o/- R o r ) 2 /(2T) 


( 13 ) 


The factor c* = 1 — i/3o/(2TR 0 ) rotates the integration path for the distorted wave functions 
ip p into the complex plane. 

Using Eq.(^) we obtain the correlation function from Eq.(||). Since the system is spher¬ 
ically symmetric we investigate the correlation as a function of the the momentum-difference 
vector ci ou t pointing in the direction of the average pair momentum K and q s ^ e being or¬ 
thogonal to K. It turned out that the correlation function C '2 can well be approximated 
by 

c 2 = 1 + exp (-qLe^Le - q L#L)> ( 14 ) 


where R S id e and Rout characterize the extension of the source in and perpendicular to the 
direction of the pair momentum K. 


IV. CASE OF ZERO CHARGE 


Before we turn to the full problem let us discuss the case of a negligible potential U. 


The distorted waves in Eq. (13) can be replaced by plane waves 


[ip p (ar)]*ip p '(otr) = exp [—iRe(a )qr + 2/m(a)Kr] 


(15) 


using the definitions ©■ To discuss the main effect it is instructive to find an analytical 
approximation. Expanding the integrand in Eq. ( |l3| ) for small flow velocities up to order 
Pq and using u — u' = Kq /to for small momenta q one obtains from Eq. (||) the result 


C 2 = 1 


exp < - q 


(14 


Rl 


1 + Pq rh/T 


+ (^) 2 t 2 

qm 


(16) 


with rh = y/m 2 + K 2 . 

Comparing the last equation with Eq.(|ld|) one identifies the first expression in the round 
bracket with the sideward radius while the whole square bracket is the outward radius with 
|K|/m being the pair velocity. The value of \a\ 2 exceeds unity by at most a few percent 
for realistic values of flow velocities. Eq. © contains the well known result |jj that the 
radial flow reduces the observed radius R S ide with increasing pair momentum and decreasing 
temperature. From the approximate behavior of the radius Rq/\J\ + P'^rh/T one recognizes 
that the dependence on the pair momentum is essentially a relativistic effect that is caused 
by the change of the relativistic pion mass to. Eq. © also contains the fact that the 
outward radius is larger for a finite pion emission time r. However, effect of mean fields 
could violate this statement as shown in the following. 
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V. NUMERICAL TREATMENT 


To incorporate the mean field we solve the Klein-Gordon equation 

o 2 

[ - 7^3 - (\/to 2 + p 2 - U ) 2 + m 2 ] (r) = 0. 


(17) 


The boundary conditions are chosen such that ip behaves asymptotically like an outgoing 
wave in the direction of p with incoming spherical waves. This is indicated by the upper 
index (-). We mention that functions with outgoing spherical waves can also be used to 
calculate the matrix element ([|) applying the relation ipp ^ = ip^p *. 

The potential U 


u = ±z e -<t( ' ll 

r W2R 0 to 2A 


(18) 


contains the Coulomb potential of the Gaussian source 0 ) with charge number Z , and 
the quantity $ denotes the error function. Further we include the possibility that the pions 
might be absorbed within the source. For this purpose we have also introduced an imaginary 
part which depends on the mean free path A. The positive imaginary part ensures that the 
wave function ip p 1 increases in the direction of the outgoing momentum p. The potential 
(|l|) is only a rudiment of the standard pion potentials used in calculating pion-nucleus 
scattering and is usually derived as part in a Schrodinger equivalent equation, see e.g. refs. 
00 - However, the simple form in Eq. ( 0 ) suffices for the study of the effect of opaqueness 
which arises from the absorption and reemission of pions in the matter. 

The solution is numerically obtained by expanding the distorted wave into partial waves 
fi 


^p } ( r ) = — J2 il(i 1<TI fl( r ) Y lm(P) Y lm(r), 

l m 


(19) 


where the quantities Y[ rn denote the spherical harmonics and the symbols 07 are the Coulomb 
scattering phases. In order to use the non-relativistic standard method known from optical 
model calculations the numerical integration was extended to large radii Rmax to render 
the term U 2 in Eq. 0 ) negligible. Once the correct radial function has been obtained 
it is analytically continued by integrating the radial differential equation from r = Rmax 
to r = ctRmax- The value obtained at this point is used to normalize the function ip(ar) 
obtained by integrating the radial equation along the path ar. 


VI. RESULTS AND DISCUSSION 

We study the model for a situation which is typical for a collision of Au on Au nuclei at 
bombarding energies of 1 GeV per nucleon. In nearly central collisions a system of charge 
Z = 120 is formed which a temperature of about T = 80 MeV for pions and a flow velocity 
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of about (/?) = 0.32 |l8| , |l9|| . We use a source radius of Rh = 10 fm which corresponds to the 
parameter i?o = Rh/'J 5. 

The correlation function (||) has been calculated for various pair momenta K as function 
of the relative momentum q. In all cases the obtained correlation function has nearly a 
Gausian shape and is fitted to Eq.(|l4|) in the region of C 2 = 1.5 to obtain the HBT radii 
Rout and Rside • 

In Fig. 1 the ratios of the fitted source radii to the true radius are displayed as a function 
of the averaged pair momentum |K| for two extreme values of the flow velocity (/3) = 0 and 
(P) = 0.5. Compared to Rq the observed radii R S ide are increased for negatively charged pion 
and diminished for positively charged pions in comparison to the true radius. An opposite 
but smaller effect is seen for the radii R ou t- These changes are significant for pions with 
momenta below 300 MeV/c. Results for pair momenta smaller than 100 MeV are not shown 
since the stationary approach is not justified. Below the Coulomb threshold the behavior 
of the radius changes drastically, see ref. |j;] for details. For comparison we have inserted 
into right hand panel of Fig. 1 the sideward radius extracted for 7 r° mesons. This curve 
nearly averages the lines for the unlike charged pions. Comparing the curves for the two 
flow velocities one recognizes that the corrections arising from the Coulomb field and the 
flow field add up nearly independently. Fig. 2 shows the ratio Rhbt/Rq for a duration of 
the emission of r = 4 frn/c which leads to the expected increase of the outward radius. 

Differences between extracted HBT radii for positively and negatively pions have been 
measured at AGS Q in the projectile rapidity region. A ratio of R S ide{^~ 7r _ ) to R S ide{ 7r + 7r + ) 
of (5.6±0.7) fm / (3.9±0.8) fm = 1.4±0.3 has been found which agrees with our predictions. 
The ratio of the outward radii of (5.8 ± 0.5) fm / (6.5 ± 0.5) fm = 0.9 ± 0.15 is smaller than 
unity although the relatively large error bars do not allow a definite comparison. These 
measurements also qualitatively agree with the ratio of 1.2±0.4 of the sideward radii observed 
by Pelte at al. 0 . In those measurements the same increase was found for the radius R ou t 
contradicting our predictions. 

Now we investigate the effect of opaqueness of the source. In ref. || a drastic change 
of the the HBT radii was predicted. Pions with momenta around ko = 270 MeV/c have a 
large total cross section with nucleons exciting strongly the A resonance. Therefore, those 
pions have their last interaction points within a thin surface zone near the direction of their 
momenta. The thickness of this zone is determined by the mean free path A of the pions. 
The inverse path length is estimated to be 

1 . _± N _± Z 

Kl = » K 7 + ", j) (20) 

which is proportional to the baryonic density n and the total cross sections crjj of pion- 
nucleon collisions averaged over the thermal motion of the nucleons. Due to isospin coupling 
the cross sections a* = are by a factor of three larger than the remaining two. This 

creates an isospin asymmetry of the mean free path in neutron rich matter. The thermal 


motion widens the A resonance to F = 240 MeV and reduces the maximum cross section 
by about a factor of two resulting in a value of 100 mb for <t£ . At normal nuclear matter 
density of no = 0.16 fm~ 3 one obtains for heavy nuclei like Au values of A° + = 1.05 fm and 
A°_ = 0.85 fm. Similar values are known from BUU calculations [Q. Now we can simulate 
the opaque source by introducing the momentum and density dependent mean free path 
A = A°no/n(l + (2 (p — fc 0 )/T) 2 ) into the potential JTs|). 

Fig. 3 shows the effect of the opaqueness for the Gaussian density distribution © 
with Rh= 10 fm for positive and negative pions without considering the electric charge. The 
opaqueness increases the radius R S ide while the radius R ou t decreases as a consequence of the 
relatively thin middle part of the half-moon shaped source region ||. The curves reflect the 
resonance shape of the absorption. An essential effect of the opaqueness is that the difference 
Rout — Rside could become negative which may compensate the positive contribution from 
the emission time. In the ultra-relativistic regime one should also add the effect of the 7r _ 
7 r + scattering since the pion density is large and the cross section could reach values up to 
15 mb for pion momenta around 200 MeV/c. 

A strong dependence of the side-correlation on the pion-pair momentum has been found 
for the collision of Au on Au at a bombarding energy of 1.06 GeV per nucleon in ref. ||. 
A comparison to that data give us a good opportunity to illustrate how the different effects 
discussed so far change the true radius. We employ the parameters T = 80 MeV and {(3) = 
0.35 as before, however we reduce the central charge to Z e ff = 60 to diminish the Coulomb 
effect trying to correct for the expansion during the pion emission. Using a hard sphere 
radius of Rh = 8 fm and an emission time of r = 4 fm/c the obtained sideward and outward 
radii are shown together with the measurements || in Fig. 4. 

We do not intend to fit the data to our parameters since our simple model lacks essential 
features, especially the time evolution of the collision treated in recent dynamical models. 
Fig. 4 shows however that essential deviations from the true radius are to be expected and 
high precission measurements are needed to gain insight in the dynamics of nuclear collisions. 
For negative pions it is clearly seen that the measured sideward radius depends stronger on 
the pair momentum than one would expect from calculations using a fixed radius. This 
means that indeed the fast pions come from an earlier more compressed stage of the matter 
while the low energy pions are emitted from a zone with a larger size. 


VII. CONCLUSIONS 

The nuclear Coulomb field increases (decreases) the observed HBT radii extracted from 
sideward correlations for negatively (positively) charged pions. The influence is opposite 
and smaller for the outward correlation. This effect is the largest for small momenta and 
is superimposed on the overall reduction caused by the radial flow. The opaqueness due to 
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pion rescattering leads to a decrease of the outward radius and an increase of the sideward 
radius. The decrease could compensate the general increase of the outward radius due to 
the duration of the emission process. 
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Figures 



Fig. 1. Ratios of sideward and outward HBT radii to the true radius Rq of a Gaussian 
shaped source as a function of half the pair momentum K. The ratios have been calculated 
without (left panel) and with (right panel) radial flow of mean velocity {(3). 
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Fig. 2. Ratios of sideward and outward HBT radii to the radius Rq of a Gaussian 
shaped source with radial flow for an emission time of 4 fm/c. 


12 








2 


1.5 



1 


0.5 


0 


pion opaqueness 

T= 80 MeV 

side correlation 

out correlation 
t = 0 fm/c 


Z = 0 
<( 3 > = 0.32 


_ 71 TV 

_ 7T + 7T + 


_ 7T TV 

_ 7T + TT + 



200 400 600 

K(MeV/c) 


Fig. 3. Ratios of sideward and outward HBT radii to the true radius Rq affected by 
flow and opaqueness of the source as a consequence of the small mean free path of pions 
within the source. 



Fig. 4. Comparison of measured sideward and outward radii || to calculated HBT 
radii affected by flow, central Coulomb field and opaqueness as a function of half the pair 
momentum. The large extension of some of the horizontal error bars indicates that the full 
range of pair momenta has been used to extract the radius. The calculation are carried 
out with a fixed source radius of Rh = 8 fm which cannot fully explain the momentum 
dependence of the measured sideward radius. 


13 
















